Spin dynamics of the doped t-J model 
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Using the exact diagonalization technique we study the low energy spin excitations in moderately 
doped (hole concentration < 25%) finite clusters of t — J model. To clarify whether a given low 
energy spin excitation corresponds to a particle-hole transition in the 'quasiparticle band' near the 
Fermi energy, we compare the electron addition and removal spectra of the respective final state 
to that of the ground state. We find that the low energy spin excitation spectrum is composed of 
different branches, with the dominant low-energy spin excitation at (n, tt) being a spin-wave like 
collective mode, whereas the low energy spin excitations with momentum transfer other than (tv, tv) 
correspond to particle-hole excitations. 
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The spin dynamics of strongly correlated electron 
systems is a key issue for the understanding of high- 
temperature superconductors. As an experimental fact, 
neutron scattering experiments show a broad peak in the 
dynamical spin correlation function for in-plane momen- 
tum transfer (jr, tt) at low frequencies ~ 40meV [jlj. Ex- 
act diagonalization studies of small clusters oft— J model 
usually reproduce this behaviour |^|,^) , and it is the pur- 
pose of the present work to distinguish between different 
interpretations which are possible for this peak: within a 
weak coupling-like picture, the dynamical spin suscepti- 
bility is determined by particle-hole excitations of either 
real electrons or 'spinons', the strong scattering intensity 
for momentum transfer (tt, tt) then would be the conse- 
quence of Fermi surface nesting. A different interpreta- 
tion would be to attribute the peak to the remnants of the 
short wavelength spin waves present in the undoped com- 
pounds; as long as their wavelength is shorter than the 
spin correlation length, such spin wave like excitations 
might well persist, but being no longer the Goldstonc 
modes associated with antiferromagnetically broken sym- 
metry, a finite frequency would be naturally expected [Q . 
Our results show that in finite t—J model clusters (lattice 
size < 20) with low hole concentration the second inter- 
pretation is the correct one: the dominant low-energy 
peak at (tt, tt) corresponds to a collective spin wave-like 
mode, particle hole excitations are restricted to different 
momentum transfer. The t—J model reads: 
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Here the Si are the electronic spin operators and the 
sum over < i,j > stands for a summation over all pairs 
of nearest neighbors. The operators Ci jCr are expressed in 
terms of ordinary fermion operators as Cj )(T (l — rii t - a ). 
To clarify the nature of the low-energy spin excitations, 
we applied the following scheme: using the Lanczos algo- 
rithm we computed the dynamical spin correlation func- 
tion (SCF) 
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where \^fu,n) {E v ,n) is the V th eigenstate (eigenenergy) 
with n holes (in particular v=0 implies the ground state). 
Fig. [I] shows S(q, u>) in the 18-site cluster with 2 and the 
16-site cluster with 4 holes (in this paper we restrict our- 
selves to these systems, whose behaviour is representa- 
tive). We then extract the energies E v , n of the dominant 
low-energy peaks indicated by arrows in Fig. 0, these 
are the two lowermost peaks for each cluster. Each of 
them is a single, unsplit peak so that a unique (and for 
the low energy peaks highly precise) value for E v>n can 
be obtained. Using the inverse iteration algorithm || 
we can now converge out the corresponding final state 
wave function l^^.n) (using a reasonably accurate esti- 
mate Etr for an eigenvalue, the inverse iteration method 
essentially consists in applying powers of (H — Etr)^ 1 to 
some trial state) . Thereby we also obtain an independent 
estimate for E^ n , which always agreed with the value 
extracted from the correlation function to an accuracy 
better than 10 -10 , essentially the limit of the Lanczos 
procedure; total spin and point group symmetry of the 
obtained |$ v ) were cross-checked to be compatible with 
the selection rules. In the last step we calculated (again 
via the Lanczos method) the electronic spectral functions 
A^(k,tu) = A^lik^uj) + A^lik,^), where 

4 Q l(fc,c) = ^|(vI/ M ,„ +1 |c fc , c7 |* Q )| 2 

5(uj — (E^.n+i — Eo >n )), 

AW(fe, W )=^|(^ )n _ 1 |4 !(T |* Q )| 2 
n 

5{lo - (£„,„_! - E , n )). (2) 

where |^ a =o) denotes the n-hole ground state and |^ r Q =i) 
the SCF-final state. These functions are simply the 
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photoemission (PES) and inverse photoemission (IPES) 
spectra for the ground state and SCF-final state. By 
inspection of the difference A^(k,u>) = A^ (k, uj) — 
A^(k,Lu) (which henceforth will be referred to as the 
'difference spectrum for the state I'iV.n)') we can now 
decide whether the creation of a spin excitation is asso- 
ciated with the transfer of single particle spectral weight 
in the band near the Fermi energy or not. 
To begin with, Fig. || shows the single particle spec- 
tral function of the ground state and the difference spec- 
trum for the SCF final state for momentum transfer 
q — (27r/3, 27r/3). In order to eliminate shifts of spec- 
tral weight on small energy scales, a relatively large 
Lorentzian broadening of OAt is used. We locate the 
Fermi energy as separating the highest PES and low- 
est IPES peaks in the ground state spectra. The dif- 
ference spectra have remarkably small weight as com- 
pared to the original spectra, particularly the parts re- 
mote from Ep: small 'wiggles', indicate slight shifts of 
'peaks' within the PES spectrum. Major changes of the 
spectral function are restricted to the neighborhood of 
E F : at fei = (-2tt/3,0), PES weight is shifted to the 
IPES side right at Ep, indicating that at this momentum 
the SCF-final state has a reduced electron occupancy of 
the quasiparticle band as compared to the ground state. 
At fci + q = (0, 27r/3) the reverse can be seen, i.e. the 
SCF-final state has an enhanced electron occupancy of 
the quasiparticle band. We thus see a shift of low en- 
ergy spectral weight with momentum transfer q, precisely 
the expected signature of a particle-hole excitation. To 
a lesser degree, such a shift can also be seen between 
(7r/3, — 7r/3) and (7r,7r/3), but the weight transfer is sub- 
stantially smaller. 

The particle-hole character is even more pronounced for 
the low energy spin excitation with momentum transfer 
q = (-7t/2, 7r/2) in the 16-site cluster with 4 holes (Fig. 
^) where we clearly see the transfer of an electron right 
at Ep from (71", 0) to (7r/2, — 7r/2). We conclude that our 
technique indeed provides a sensitive tool for detecting 
particle hole excitations and consequently apply it to the 
dominant low energy peak at (ir, tt). The difference spec- 
tra for the respective final states in the 18 and 16-site 
clusters are shown in Figs. [| and [|. There are 'wiggles' 
near the Fermi energies in some of them which indicate 
a shift of the quasiparticle peaks. There is, however, no 
appreciable shift of weight from PES to IPES near the 
Fermi energy, despite the fact that the spectral intensity 
of these excitations in S(q, w) is much stronger than the 
particle-hole excitations considered previously. The spin 
excitation at (it, tt) thus obviousy is not associated with 
a shift of spectral weight near the Fermi energy, i.e. it 
is not a particle-hole excitation but rather a collective 
mode the character of which remains to be clarified. 
Up to a constant, the spin correlation function for this 
momentum equals the correlation function of the stag- 
gered magnetization, Ms = J2 ie A S i ~ J2ieB S i- Wc 
will now assume that the ground state expectation value 
(Mj) = (N ■ ms) 2 , where N denotes the number of sites 



and ms is of order unity, and show that this allows the 
explicit construction of a low energy collective spin ex- 
citation with momentum transfer (tt, 7t) (we do not as- 
sume (Ms) 7^ 0, i.e. there need not be broken sym- 
metry). We decompose the ground state wave function 
l^o) as I'I'o) = X^S^lVv), where the are defined 
by the requirements Ms|i/>„) = v\^) v ) and ($> v \tyv) = f- 
In other words, we decompose the ground state into 
components with fixed value of the staggered magneti- 
zation. We define the matrix elements of the Hamilto- 
nian, h u ,u — (ip v \H\tpu), and obviously have h V n = 
for \u — fi\ > 2 (the Hamiltonian can change the stag- 
gered magnetization at most by 2). The requirement 
that the ground state energy Eq be extensive then im- 
plies h^u ~ N. 

Introducing the translation operator by one lattice site, 
T10, the translational invariance of \^q} guarantees that 
5i/Ti,o|V , f) = 9— u\">P-u)- We now define the trial state 
= Z)w>o(^f IV'w) - 9-v\ip-v))- By construction this 
state changes sign under T^o, i-e. it has momentum 
(tt,tt). The expectation values of the Hamiltonian with 
|* ) and differ by AH = h , g 2 + 4fti,o5i5o + 

4/i2,o52<?o + +4ft-i,-iffi5-i, the norms of these state differ 
by AN = g 2 . The energy A = (AH + E a AN)/(l - AN) 
thus provides a rigorous upper bound for the energy dif- 
ference between ground state and lowest state with total 
momentum (tt, 7r) and restricting the Hilbert space to 
|*o) and we find S(Q, w) ~ m 2 s S(uj - A). 
Fig. pi shows different forms of the probability distribu- 
tion g^ which are compatible with the requirement that 
( M s) = Yjv v2 9 2 v ~ n2 - Form ( a ) implies that g v = 
for v ~ so that A=0. Hence |\&o) and are degen- 
erate and e.g. the state |\I>+) = 5Z J/>0 S^IVv) by itself 
is already an eigenstate: this corresponds to true bro- 
ken symmetry. Form (b) implies a nonvanishing A; here 
the number of g 2 which differ from must be of order 
N, so that g v ~ 1/VN and hence A ~ N°: we have a 
potentially low energetic collective mode in the spin ex- 
citation spectrum for all system sizes, which corresponds 
to a change of the relative phase of the two components of 
the ground state wave function with positive and negative 
staggered magnetization. If g 2 continuously approaches 
the broken symmetry form (a), the excitation energy of 
this mode will approach zero. For completeness we note 
that (Mj) ~ N (as is the case for a free electron gas) 
would imply A ~ yN, so that we cannot get a meaning- 
ful low energy mode from the above construction. 
Let us now check this interpretation of the (it, it) 
mode. The exact diagonalization technique gives 
I'I'o) directly as a linear combination of basis states 
cj „ cj „ „ \vac), so that it is easy to construct 

the state \^\) (normalized to unity). We can then com- 
pute its overlap with the SCF final state for momentum 
transfer (tt, it), \%>scf), i-e- the quantity | (*i |*scf) | 2 - 
Results for various cluster sizes and dopings are summa- 
rized in Tab. |, which quite obviously confirms our hy- 
pothesis: there is always a substantial overlap between 
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the respective |<J/i) and the exact SCF final state, and for 
hole concentrations around 10% (2 holes) this overlap is 
quite comparable to that in the undopcd Heisenberg an- 
tiferromagnet, where we can expect that the (tt, it) mode 
will evolve into the spin wave mode for infinite system 
size. The finite energy collective mode, which we could 
predict on quite general assumptions thus indeed seems 
to be realized in the small clusters. From our general 
contraction we expect that this mode persists also in the 
infinite system and thus may explain the neutron scat- 
tering experiments. 

In summary, our results show a fairly complicated pic- 
ture of the spin excitation spectrum in the doped two- 
dimensional t— J model. For low doping levels, (i.e. 2 
or 4 holes) the dominant low energy mode is located at 
(it, it), and we have identified this mode as a collective 
mode comparable to the spin waves present in the un- 
doped cluster. Away from (tt, tt) on the other hand, we 
have fairly conventional particle- hole excitations. 
The fact that parts of the spin excitation spectrum of the 
Heisenberg antifcrromagnet persist in the doped system 
appears very natural if one adopts the picture that the 
doped system should be modelled as an incoherent 'spin 
background' in which spin-bag-like quasiparticles corre- 
sponding to the doped holes are moving ||. Then, it 
seems natural to distinguish between two types of spin 
excitations: the first are the particle-hole excitations of 
the hole-liquid, which should resemble that of a Fermi 
liquid with a Fermi surface volume proportional to the 
number of holes. This part of the excitation spectrum 
can quantitatively explain the Pauli susceptibility at low 
temperatures j7|]. In addition there are the excitations 
of the 'spin background' which may, with some modifica- 
tions, resemble that of the Heisenberg antiferromagnet. 
One example here is the spin wave-like mode at (tt, tt), a 
possible other example the 'two magnon' excitation ob- 
served in the Raman spectrum until well in the super- 
conducting doping regime || . 

To conclude, we note that the general scheme developed 
above allows for a highly selective study of the character 
of arbitrary low energy excitations; it thus seems possible 
(at least in principle) to obain a fairly complete picture 
of the entire low-energy excitation spectrum of the 2-d 
t — J model clusters. 
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FIG. 1. Spin correlation function in the 18-site cluster with 
2 holes (top) and the 16 site cluster with 4 holes (bottom). 



FIG. 2. Difference spectrum for the SCF final state in the 
18-site cluster with 2 holes with momentum (4^, ^). The 
black (gray) area corresponds to IPES (PES). Also shown is 
the single particle spectral function for the ground state (full 
line: PES, dotted line: IPES). The vertical line marks the 
Fermi energy. 



FIG. 3. Difference spectrum for the SCF final state in the 
16-site cluster with 4 holes with momentum (~, £). All con- 
ventions are like in Fig. ^| 



FIG. 4. Difference spectrum for the SCF final state in the 
18-site cluster with 2 holes with momentum (n,n). All con- 
ventions are like in Fig. ^| 



FIG. 5. Difference spectrum for the SCF final state in the 
16-site cluster with 4 holes with momentum (n,ir). All con- 
ventions are like in Fig. ^| 



FIG. 6. Two possible forms of the 'wave function' g v . (a) 
implies broken symmetry, (b) implies no broken symmetry. 



TABLE I. Overlap of the exact SCF-final state and the 
trial state based on the spin-wave picture for different cluster 
sizes and dopings. 20 sites and 4 holes exceeds our computer 
capacity. 





16-site 


18-site 


20-site 


holes 


0.8378 


0.7560 


0.8214 


2 holes 


0.7187 


0.7388 


0.7059 


4 holes 


0.4346 


0.5620 
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